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ABSTRACT 


This thesis examines parametric x-radiation (PXR) which is the Bragg scattering of the 
virtual photons associated with the Coulomb field of relativistic charged particle from the 
atomic planes of a crystal. Higher order parametric x-radiation from the {002} planes of a 
thick mosaic graphite crystal have been observed. The raw PXR data was collected using 
a SiLi detector and a Pulse Height Analyzer (PHA) software program. The data was 
corrected for various effects including attenuation, detector dnft and efficiency. The 
absolute number of photons per electron was obtained by using the fluorescent x-ray yield 
from a tin foil backing on the graphite crystal to determine the LINAC current. The 
number of photons per electron observed greatly exceeds the expected values. 
Comparison of the ratio of intensity of a given order to the first order [I(n)/I(1)] to the 
theoretical ratio shows that the ratios increase with order. Not only is the absolute 
intensity greater than expected, but the higher orders (compared to the first order) are 
larger than expected. Lastly, the intensity for various crystal angle orientations and a fixed 


detector angle was measured. 
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I. INTRODUCTION 


Parametric X-radiation (PXR) 1s the scattering of the virtual photons associated with 
the Coulomb field of a relativistic charged particle by the atomic planes of a crystal at or 
near the Bragg condition. It is monochromatic and spectrally intense. Potentially, PXR 
may be an important new X-ray source for a host of applications [Ref. 1]. Ter-Mikaelian 
(1972) first developed the theory of PXR as resonant radiation produced in a thin crystal 
(Ref. 2]. Later Feranchuck and Ivashin developed a theory for thick crystals. In 1985, the 
first experiment verifying this theory was performed by Baryshevsky et al at the Tomsk 
synchrotron [Ref. 3]. This experiment used a 900 MeV electron beam, bombarding a 
diamond crystal at a Bragg angle of 45 degrees producing x-rays in the 5 to 25 keV range. 

The first experiments to observe PXR outside the USSR were conducted at the Naval 
Postgraduate School (NPS). [Ref. 2]. The first experiments were done with 20um Si 
(Ref. 2], 44, 320 um Si [Ref. 4] and Carbon [Ref. 5]. 

This thesis is an in depth study of the enhancement of the higher order PXR due to the 
production and subsequent attenuation of x-rays in a thick mosaic carbon graphite crystal. 
This work presents measurements of absolute PXR production efficiency and the effects 
of experimental apertures on the bandwidth of PXR spectral line widths, and makes 
comparisons with calculations. Effects due to beam crystal alignment and due to limited 
detector aperture field of view have also been investigated. These crystals reveal greater 
relative production rates for the higher order x-rays than can be explained by target 
thickness effects. Furthermore, these measurements show that PXR production for all 
orders in a thick mosaic graphite crystal greatly exceed predictions from calculations based 


upon current theoretical understanding. 


Il. THEORETICAL BACKGROUND 


Parametric x-radiation (PXR) is one of several types of x-radiation that are produced 
by the interaction of a relativistic charged particle propagating in structured media. These 
include bremsstrahlung, transition radiation (TR), channeling radiation (CR) and 
Cherenkov radiation. These various mechanisms are discussed and compared in reference 
6, which claims that PXR is the most effective source of quasi-monochromatic and 


focused x-rays, and reference 7 which compares PXR to synchrotron radiation. 


A. PARAMETRIC X-RAYS (PXR) 

Parametric x-rays are produced when the electromagnetic field of the relativistic 
charged particle interacts with the perodic electric susceptibility associated with the 
crystal lattice structure. When the charged particle beam satisfies the Bragg condition 
with the crystal lattice, virtual photons are diffracted as real x-rays. [Ref. 2] 

PXR is produced when the photon index of refraction n(K,@) is greater than unity 
when the Bragg condition 

ke =(K, +7)’, (1 =|x,|) (1) 
is fulfilled, where Tt is the reciprocal lattice vector, kK is the wave vector of the emitted 
photon, and @ Is its frequency. This condition is illustrated in Figure 1. As a result, the 
condition for Vavilov-Cherenkov radiation 

1- =" n(x) =0, (n=c=1) (2) 
is satisfied in a crystal. This is in contrast to a homogeneous medium, where the index of 
refraction is less than unity for x-rays, and Vavilov-Cherenkov radiation is not produced. 
Therefore, PXR is simply Vavilov-Cherenkov radiation produced by diffraction according 


to equation (1). [Ref. 8] 
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Figure |. Diffraction of virtual photons associated with the charged particle beam by the 
atomic planes of the crystal showing the Bragg angle as the angle between the charged 
particle beam and the crystal planes. [Ref. 2] 


There are two ways to treat the diffraction of real and virtual x-rays: 1) dynamic theory 
and 2) kinematic theory. The first case is realized in ideal crystals of thickness L, greater 
than the X-ray extinction length. On the other hand, the kinematic theory of PXR is 
applicable for real crystals, which consist of thin mosaic blocks turned relative to each 
other at an angle d>nV/E, where m is the mass of the charged particle and E is its energy. 
The kinematic diffraction is more convenient for the experimental observation of PXR and 
it can be most easily analyzed because the angular and spectral distributions of PXR are 
essentially simplified and have a universal form for different crystals. [Ref. 8] 

Ter-Mikaelian formulated the expression for the number of photons emitted by a 
charged particle in a thin crystal; treating the interaction as a perturbation, with the crystal 
requirement 

KL|n -1| <1, (3) 
where L is the effective thickness of the crystal. If the refraction of photons and the 


multiple scattering of the electrons in the crystal are taken into account, then the results, 


for a mosaic crystal with 6>m/E and with PXR from different blocks being coherent, 


remain correct. The spectral and angular distributions of emitted photons are formulated 
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Here gr and gg are the Fourier components of the dielectric susceptibility parallel and 


perpendicular to the reciprocal lattice vector. These values are directly related to the 
coherent scattering amplitude of the photons by the atoms of the crystal. 6, is the angle of 
multiple scattering. Equation (4) 1s wnitten in a coordinate system with the z-axis directed 


along the particle velocity ¥; K, =an,; and L, = [wIm(g,)] (Ref. 8] The differential 


radiation cross section for PXR as derived by Feranchuck and Ivashin can also be wnitten 


in terms of angular displacements from the Bragg angle instead of n, as: 
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x-ray scattering factor of the crystal, V is the volume of the unit cell of the crystal, re is 


the classical electron radius, fi=c=l, p, =A/ (L,(w)22) and L,(w) is the photon 


absorption length. [Ref. 8] 

The characteristic quasi monochromatic photons produced at large angles relative to 
the charged particle's velocity vector, v, distinguishes PXR from other radiation 
mechanisms. The particle energy does not determine the emitted photon's wave vector, K, 
nor frequency ®@p. The emitted photon’s wave vector depends upon the reciprocal lattice | 
vector, T, the velocity vector of the electron and the Bragg frequency. The Bragg angle, 
Op, is the angle between the velocity vector, v and the crystallographic planes associated 
with the reciprocal lattice vector t. The PXR maxima are defined by the crystal lattice 
structure amplitudes corresponding to tT. The values of the emitted wave vectors, kK, and 


frequencies Wp are 


K=0,V +7; oN = Pike (7) 


a sais sees)... (8) 





where d is the distance between the crystal planes corresponding to the vector T. [Ref. 3.] 

The intensity of PXR is sapedtional to lg} and is greater in crystals with higher 
packing factors, such as a diamond lattice. Rotation of a crystal by some A@ results in the 
emitted wave vector, kK, rotating 2A0 in the same direction while the frequency varies as a 
result of equation (7). The angular distribution of emitted photons can be derived by 


integrating equations (4) or (6) over w, resulting in 
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where m is the rest mass of the electron, Lg 1s the absorption coefficient which is a 
function of @ and 8poh is the characteristic angular spread of the PXR. 9, includes the 
effects of multiple scattering of the beam electrons, electron beam divergence and crystal 
mosaicity. This can be put in dimensionless form by means of the normalized amplitude J= 


N/No and angles x,y = 9x y/Sph; 
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The number of photons detected in an angular width of 8, about 9, is obtained by 


integrating equation 9: 


N, = mN,(1+cos° 204) Pap. (14) 
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where Pp, =—*. The detector angular size is a factor in determining the value of Np, 


ph 
even for 8p<<®8pp. This factor does not occur in channeling or in bremsstrahlung 
radiation; however in PXR, the intensity slowly decreases inversely proportional to 9 as 8 


increases. [Ref. 8] 


According to equations (14) and (15), there is no threshold characteristics in the PXR 


intensity. The angular frequency, however, depends on Oph(E) which increases with the 


electron energy 


Eben (16) 





in accordance with equation (9) and (10), thus setting the limiting threshold energy for 
PXR. The number of photons decreases rapidly for E<E opt. The number detected 


becomes 


opt 
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This phenomena suggests that if the beam energy 1s not sufficient, regardless of whether 
all the other parameter conditions are met, PXR may not be observed. [Ref. 8]. 


Integration of equations (4) or (6) over n, or dQ respectively results in the frequency 


distribution; 
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Equation (18) can be put in terms of frequency by using the chain rule a = gia and 
OW du 0M 


= can be determined using equation (19). Although equations (18) and (19) do not 


apply to real experiments with finite aperture, they have been previously applied to the 


experimental results. [Ref. 8] 


B. MULTIPLE SCATTERING, BEAM DIVERGENCE AND MOSAICITY 


In reference 1, the Adischev et. al. have used 


FE Pae 
C=. = — | 20 
Ms SET (20) 
E,=21 MeV, L is the effective thickness of the crystal and Lp is the radiation length. This 
value is an approximation applicable to very high energy electrons and is well known to 


Over estimate the multiple scattering effect. [Ref. 9] The Bethe-Ashkin formula for 


multiple scattering which gives a more accurate estimate for 6;,, is used in the analysis 


and is given by; 
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where Z is the atomic number, A is the atomic weight and t is the thickness of the foil in 





(21) 


grams per square centimeter. The charge of the scattered particle is z, p 1s its momentum 
and v its velocity, <@2> is an average scattering angle value obtained by integrating over 
angle, where 6) 1s the maximum angle used in the analysis. It is estimated as the angle for 
which there 1s, on the average throughout the full thickness t of the foil, only one collision 
with 8 > 0). [Ref. 11] 

Since x-rays are produced throughout the electron's path length, the amount of multiple 
scattering varies from zero to the value appropriate for the full path length of the electron 
in the crystal. Additionally, the x-rays which are generated inside the crystal are 
attenuated due to absorption. Thus there is a weighting of x-rays reaching the detector, 
with the weighting factor equal to e", where yu is the absorption coefficient for the 
particular x-ray energy. Therefore, x-rays produced near the front face of the crystal are 
weighted more than those near the back surface. This weighting effectively biases the 
detected x-rays toward those produced at lower values of x, i.e. near the front face of the 


crystal where the multiple scattering is lowest. The proper way to account for the effects 


of multiple scattering would be to add up all the x-rays produced as a function of x, 
including an x dependent multiple scattering effect. 

For simplicity, a mean value of <®,4 5 2> can be calculated by weighting the value of 
<@“> with the absorption attenuation of x-rays, e"”. The mean value for x is defined as 
follows 

xe dx 


< x >= —_—__.. (22) 
| Cia 


0 
[Ref. 9] Since the Bethe-Ashkin formula for <®y4 s 2> to first order is proportional to x, 
this value for x is used in equation (21) in place of the crystal thickness, t, to determine the 
multiple scattering angle 8, [Ref. 11]. In our treatment, the effect of the divergence of 
the electron beam is not included. To more migorously include the effect of beam 
divergence on the x-ray angular distribution and bandwidth, the electron beam angular 
distribution must be convolved with the angular distribution of the x-ray production as has 
been done e.g. in the case of transition radiation [Ref. 12]. For simplicity of calculation 


the effect of divergence was not included in the analysis since the effective beam 


divergence of the NPS LINAC (<1 mrad) is unimportant. 


C. EFFECTS OF FINITE APERTURE 

Previous PXR experiments to date have compared the experimental data with the 
results of equation (18). However, the theoretical values obtained by using equation (18) 
refer to an unrealistic idealization i.e., an infinite detector aperture, which does not 
approximate the real conditions of these experiments and neglect the effect of the aperture 
on the PXR spectrum. 

In order to calculate the spectral distribution for real experimental conditions, equation 


(15) should be integrated over the appropriate angular field of view as determined by the 


detector size and the distance from the crystal to the detector. The integration limits 6. 


and 6, are set by the angular fields of view in and out of the plane of observation 
(A@, = 26, and A@, = 26, respectively). [Ref. 13] 
Case I. Infinite 6., finite 6. Integrating equation (7) first over Oy with a finite upper 


limit of 6, and then integrating 6, from 0 to © results in 
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where the function J> 1s defined as, 
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dependence, see equation 6. [Ref. 9] 


Case ll. 6, and 6, both finite. In this case there are two regimes. a) The first is the 


narrow line limit regime in which p, = A@_, then 


(2) 
ON a =: (25) 


b) The second is the aperture limited bandwidth regime for which the delta function 


approximation 1s still correct if A@, >> p,. Then: 
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In this case the energy bandwidth is determined based on the relative value of A@, and 
the full width at half maximum (FWHM) of J, (u). If A@, is larger than the FWHM of 
J,(u), then the bandwidth is determined to be the FWHM of J,(u). This region is 
referred to later in this thesis as the "near region". If A@. is smaller than the FWHM of 
J(u), then the bandwidth is proportional to A@, and 1s relatively insensitive to the form 


of J(u). This regime is called the small aperture regime referred to later in this thesis as 


the "far region". The bandwidth for the small aperture regime 1s: 





oo (27) 


It is the small aperture regime that applies to this thesis work and that is most applicable to 


previous PXR experiments [Ref. 10]. Figures 2 and 3 illustrate the function J, (u) with 


A@.<FWHM of J,(u)for the far region and A@,>FWHM of J,(u) for the near region 


respectively. 





Figure 2. J> and the field of view for the "far region" at the Bragg angle. 





Figure 3. J> and the field of view for the "near region" at the Bragg angle. 


Iii. PXR EXPERIMENTAL DESCRIPTION 


Several PXR experiments were conducted using a 1.39+0.01 mm thick graphite crystal 
in Bragg geometry. The x-ray observation angle was set at 45 degrees with respect to the 
electron beam direction. In this configuration, corresponding to a Bragg angle of 22.5 
degrees with respect to the {002} atomic planes of the carbon mosaic crystal, 
measurements were made at detectors placed at 100 cm. and 29 cm. from the target to 
observe the effects of changing the solid angle of observation. At both these 
configurations, "rocking data" measurements were made. "Rocking data" were obtained 
by rotating the crystal about the nominal Bragg angle of 22.5 degrees. The rocking data 
provided the dependence of the absolute and relative intensities of the various order x-rays 


when the Bragg condition was not exactly satisfied. 


A. ALIGNMENT AND ELECTRON BEAM STEERING 

The NPS Linear Accelerator (LINAC) was used to conduct these experiments. The 
beam energy varied with each experiment with a nominal electron energy of 90 MeV. 
Two different SiLi detectors [Ref. 14 & 15] were used in these experiments. The first was 
an EG&G ORTEC model 7113-16250 with an active detector region of 200 mm2. and 
the other was a Canberra model $1I200250 with an active detector region of 200 mmZ2. 
The SiLi detectors have nominal time resolution times of 12 usec determined by the pulse 
shaping preamplifier. Because the LINAC is pulsed at 60 Hz with a macro structure 
length of 1 usec, the system 1s count rate limited. The LINAC duty cycle makes the count 
rate 60 MHz when a single photon is detected during each LINAC pulse. The LINAC 
was adjusted to limit the macroscopic average count rate including background to one 
count per three to five machine macropulses in the whole detector energy range. 


Consequently, the LINAC beam was limited to a dark current estimated to be of the order 
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of 2x10-!3 amps. Dark current refers to using the LINAC with zero gun grid voltage, so 
that only stray electrons are accelerated. The ratio of counts per machine pulse was 
determined by use of fluorescent foils to prevent the double counting of photons by the 
detector. That 1s, if two photons are stopped in the detector simultaneously (within the 
peak shaping time of the detector) these two photons will appear as a single photon with 
an energy equivalent to the sum of the two coincident photons. Pulse pile up was of 
critical concern in this experiment because two first order PXR photons observed 
simultaneously would register as a second order PXR photons at double the energy. This 
would falsely enhance the higher orders as well as reduce the first order, resulting in ratios 
of intensity of higher order peaks to first order peaks too high. Maintaining a constant 
dark current was difficult to achieve and required constant operator attention. 

The experimental setup is shown in Figure 4. The beam is shown entering a vacuum 
chamber from the left at the entrance port. At the center of the chamber is a target ladder 
which can be remotely rotated, raised, lowered and tilted. An alignment-reference laser 1s 
located outside the port located 90 degrees counter clockwise from the entrance port. 
The detector is located 45 degrees clockwise from the beam exit port or 135 degrees 
counter clockwise from the entrance port. To help align the electron beam a 
phosphorescent screen was placed at the end of a 27.5 in. pipe located at the beam exit 
port. 

Experiments prior to July 22 1992 were conducted using an EG&G ORTEC Si(L1) x- 
ray detector model number 7113-16250. Due to the fixed height of the target chamber 
and the detector's large liquid nitrogen tank located below the detector window the 
detector couldn't be placed closer than about two feet from the chamber. For experiments 
with this detector a copper vacuum pipe 27.5 in. long was placed at this port fitted with a 
one mil Kapton window (1 in. diameter). The results obtained from these experiments 


will be referred to as the far case or the far region. Experiments on or after July 2294 


target chamber 





incoming 
electron 
beam 


entrance fp 





detector 


alignment laser 
Figure 4. PXR experimental setup in the LINAC end station, showing incoming charged 
particle beam and beam line, alignment-reference laser, detector, and entrance and exit 
ports. |. Target ladder with calibration foil, mirror, and carbon crystal. 2. downstream 
phosphorescent screen. 


were conducted using a Canberra Si(Li) x-ray detector model number S1200250. This 
detector consist of a two fill and vent port dewar to allow operation at any orientation and 
a "slim-line" detector chamber which is small in size and weight. [Ref. 14] This detector 
has a much smaller honzontally mounted liquid nitrogen tank and could be positioned 
directly next to the vacuum chamber. The vacuum chamber port was fitted with the same 
one mil Kapton window that covered the copper pipe. The results obtained from these 
experiments will be referred to as the near case or near region. 

A remotely vertically positionable ladder and 2 axis goniometer were used to hold and 
orient the mosaic carbon crystal used in this experiment. The crystal used in all the 
experiments 1s a mosaic graphite monochromator measured to be 1.39+0.01 mm. thick. 
manufactured by Union Carbide. The density of the crystal is 2.260 + 0.005 grams per 
cubic centimeter. The spacing of the reflecting {002} planes is 3.356 + 0.003 angstroms. 
The mosaic spread is 20.4 degrees. [Ref. 16] (Pure crystals (mosaic spread = 0) such as 
Si have lower reflectivity than mosaic crystals in normal x-ray diffraction and therefore, a 


higher the mosaic spread results in a theoretically higher measured intensity.) On June St 
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half of the crystal was covered with a 0.0275 mm thick layer of tin 7.31 g/cm on the back 
face of the crystal. This foil was essential as a calibration source. From the location of the 
Kg fluorescence line, the energy stability and energy calibration of the detector through 
out the experiment could be monitored. The observed yield of photons is proportional to 
the integrated electron current through the foil. This information was used to determine 
the absolute photon production rate. [Ref. 17] 

The ladder which housed the mosaic graphite crystal also contained three separate foils, 
copper (Cu), yttrium (Y) and indium (In), which were replaced on June sth by a sandwich 
foil consisting of a titantum (T1) foil 0.0164 mm. thick, a yttrium (Y) foil 0.0603 mm. 
thick, and a tin (Sn) foil 0.0275 mm. thick. These foils were used for energy calibration 
purposes and to determine the LINAC beam current. (The details of this calculation are 
presented in reference 17.) The ladder also contained a phosphorescent screen with a pin 
hole in the center which was used to position the electron beam. The target ladder also 
contained a mirror, coplanar with the crystal face which was used to reflect a reference 
laser back onto itself to establish the zero degrees or "home" position. 

The reference laser was used to establish the onentation and position of the crystal 
with respect to the electron beam. Using cutouts with a pinhole in the center to cover the 
ports, the laser was adjusted so that 1s went through the geometric center of the chamber. 
The ladder was then placed in the center of the chamber with the motor controls attached 
and the pinhole set to the position of the laser beam. The ladder was manually adjusted 
for tilt so that the laser went through the pinhole in the center of the phosphorescent 
screen, and remained on the pinhole upon 360 degree rotation of the crystal. A mirror and 
the graphite were positioned so that their faces were coplanar and such that the rotation 
axis of the ladder was in their planes as well. Then, by setting the mirror at 45 degrees 
with respect to the laser beam, the laser beam was reflected down the exit pipe. A second 
phosphorescent screen was placed at the end of the exit pipe with the center of the screen 


location based on the location of the reflected laser. This second screen greatly improved 
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the alignment of the experiment. It was determined that the position of the center of the 
screen was not identically at the geometric center of the pipe. After the first data were 
obtained, using this screen, the screen slipped, and could not be replace exactly without 
breaking vacuum and re-configuring it. The screen was replaced with its center near the 
geometric center of the pipe, which introduced a slight additional uncertainty in the crystal 
angle of less than two tenths of a degree. For this reason, in future experiments, OTR 
should be used on a down stream foil as an alignment as was first done with Si. [Ref. 4] 

After the chamber was sealed and a vacuum was established the LINAC was then 
tuned to the given energy, and a beam was established. The beam was adjusted so that it 
went through the center of the quadrupole magnets, so that adjustments in klystron 
frequency and phase during the experiment would not affect the location of the beam. 
Using end station steering magnets, the beam was steered so that it went through the pin 
hole located on the screen on the ladder. Since the only means of steering without 
changing the position of the beam through the quadrupole magnets is end station steering 
(see Figure 5) a compromise had to be made to get the beam through both the center of 
the screen on the ladder and through the center of the downstream screen. Additionally, 
the beams focus required a compromise. Focusing downstream leads to a more parallel 
beam at the target, but it also makes a larger beam spot. Focusing on the target gives a 
smaller beam spot, but the beam electrons are not parallel. The beam was focused 
approximately half way between the two screens based on these two factors. This 
uncertainty however, contributes to the uncertainty in the crystal angle with respect to the 
beam line of less than one tenth of a degree. 

Once the beam was focused and the steering adjusted, the electron gun was turned off 
and the current was reduced to dark current. The ladder's angle was adjusted using the 
laser since the electron beam's position is known relative to the laser. Refer to Figure 6. 

The angle 8., for positioning the crystal is determined with respect to the reference 


laser. 0, is the angle by which the crystal was rotated and was equal to 6p or Op + 180 
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Figure 5. Naval Postgraduate School 100 MeV Linear Accelerator (LINAC) showing 
klystrons, quadrupole magnets and end station steering magnets. 
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Figure 6. The relationship between the incoming particle beam and the ladder, 8p, and 
between the beam and the detector Op which ts fixed at 45 degrees. 


degrees. The home or zero degrees position was established by using a mirror on the 
ladder. When the ladder is oriented at zero degrees with respect to the laser, the beam 1s 
reflected back onto the laser. Using a remote camera the "home" or zero position is 
established and fine adjusted after every rotation of the crystal. This gives a more 


consistent angle and also removes the backlash inherent in the motor controls. With 


reference to the laser the electron beam enters the chamber at 90 degrees. The target 
crystal orientation for these experiments varied with respect to the beam line while the 


detector angle was fixed at 45 degrees. 


B. ELECTRONICS FOR SPECTROSCOPY 

Because the end station is relatively exposed to the klystron gallery (see Figure 5) the 
SiLi detectors are susceptible to ground loops and radiated noise from the klystrons. 
Consequently great effort was taken to insure sufficient grounding and shielding to 
suppress the noise pulse. In addition a gating system was set up. Figure 7 1s a block 
diagram of the electronics setup. The signal 1s received by either SiLi detector. It is 
amplified by an ORTEC Amp 571 pre-amplifier and an ORTEC 450 amplifier which was 
replaced on May 88 by a TENNELEC TC 244 amplifier. It is passed though an ORTEC 
426 linear gate with a gate width of 5 microseconds. This gate was replaced in later 
experiments with a TENNELEC TC 304 linear gate which has wider gating capabilities. 
The gate is triggered by a Stanford Research System Inc. model DG535 four channel 
digital delay/pulse generator. The detector and preamplifier were located in the LINAC 
end station. To allow for possible changes in settings the remaining components were 
located in the control room. This resulted initially in a noise problem that was eliminated 
by trial and error, using grounding straps and a return signal from the amplifier in the 
control room to an oscilloscope in the end station. The delay generator delay time (Tq) 1s 
adjusted with respect to To, the machine's start sequence time. The delay time is adjusted 
so that the start of the gate coincides approximately with the apparent arrival of the beam 
pulse at the ladder in the chamber. The ideal delay time, Tg for the gating was determined 
by using an oscilloscope along with the pulse height analyzer and the fluorescent signal 
from a copper foil. Tg was determined to be the time at which the maximum signal was 


observed. Shorter and longer delay times gave no signal. The maximum was found to be 
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Figure 7. Electronic setup. 


6 us, for the earlier experiments using the ORTEC amplifier, and 30.6 us for the 
TENNELEC amplifier. 

The LINAC generates 60 pulses of one microsecond duration every second. Since x 
rays are only produced during the actual pulse time this gating doesn't affect the PXR 
spectrum or the sandwich foil calibration spectrum, other than by limiting the noise that is 
non-coincident with PXR production. The gating limits detection to the on-time of the 
LINAC pulse. However klystron RF noise adds a bias to the pulse signal energy 
calibration which tends to be non-random with respect to Tg hence shifting the energy 
calibration compared to that obtained using radioactive sources. Since the PHA 
determines the energy (or channel) based on the energy deposited in the active detector 
region at a given time and if the x ray arrives during a consistently negative klystron noise 
the PHA determines the energy of the x ray to be lower than it is. The klystron noise was 
typically in the negative voltage portion during the gating (see Figure 8). The negative 
voltage adds to the positive voltage produced by the detected x-rays. This lowers the 
output voltage. The net effect is a lower pulse height, which is interpreted by the pulse 
height analyzer as a lower energy photon. Therefore the calibration must be obtained 
during the same portion of the machine cycle as the PXR spectrum to get the same 


klystron noise effect for both PXR and calibration x rays. For radioactive sources gating 
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Figure 8. Typical shape of klystron noise during PXR data collection. The arrow 
indicates approximate gating portion. 


effectively eliminates the signal. Without gating the noise contributes randomly and 
therefore results in a net zero shift. Gated fluorescent source's peaks and the gated PXR 
spectral peaks will consistently shift to lower values than expected if the calibration is 
based on radioactive sources. Therefore, the known fluorescent energy values were used 


to calibrate the detectors. 


C. ENERGY CALIBRATION 

For each experiment an energy calibration spectrum was obtained. These were obtained 
by measuring the x-ray fluorescence from either three separate foils or the sandwich foil 
described in Section A. Observed x-rays are the fluorescent x-rays at the K-edge of the 
excited atom. The K-edge x-ray energy for titanium 1s 4.509 keV, 8.04 keV for copper, 
14.932 keV for yttrium, 24.10 keV for indium and 25.196 keV for tin [Ref. 18] The 
calibration is obtained by establishing the centroid of the fluorescent peak. That centroid 
channel is then set equal to the known de-excitation energy. Assuming a linear 
relationship between channel position and energy, the equation for determining the energy 


for a given channel is determined by the following equations: 


E=a(chn)+8 , (28) 

eee (29) 
cn, — clin 

B=E, -—a(chn,), (30) 


where E; and E> are the energies of peak 1 and peak 2 respectively, and chn, and chn> 
are the corresponding centroid channel numbers, @ is the slope and f is the intercept of 


the Energy vs. Channel line. 
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IV. DATA AND ANALYSIS 


A. SINGLE ANGLE EXPERIMENTS (FAR REGION, A@, <FWHM OF J>) 

To obtain a PXR spectrum, the target ladder was lowered until the mirror was in line 
with the electron beam, and the home position was determined as described in Chapter III. 
section A. The ladder was lowered to the crystal and then rotated to the desired angle. 
All single angle experiments used the ORTEC SiLi detector and were in the far region. 
For the far region the detector was approximately 100 cm from the target and subtended 
2.0x10-4 sr. A summary of the various variables for these experiments can be found in 


Table 1, and the details are given below. 


TABLE 1. SUMMARY OF VARIABLES FOR THE SINGLE ANGLE EXPERIMENTS 
WITH NOMINAL ANGLES NEAR THE BRAGG ANGLE OF 22.5 DEGREES. 


Gate Detector 


a (width) | (distance to 


(peaking time) 


5 May Ortec 450 SiLi model 
(12 Us) : #7113-16250 
1006 mm 
TENNELEC i 
22. 
(12 Us) 


(1002 mm 
TENNELEC | Ortec 426 =~ 
TC 244 (35 ps) (1003 mm) 


(8 [s) [7.7 Us] 


8 July | 84.5 





"Beam energy” is the energy of the LINAC electron beam, "Amplifier" is the model of the 
amplifier, with the amplifier peaking time listed below in parentheses, "Gate" is the model 


of the linear gate used with the gate width listed in parentheses and the delay time from 


to 
to 


the LINAC start for the start of the gate are listed in brackets, "Detector" is the model 
number of the detector with the distance to the detector listed in parentheses (the 
variations are due to slight variations in the distance from the kapton window to the 
detector window referred to as the air gap), "Nominal Angle" is the crystal angle with 
respect to the electron beam line. 

Data was obtained on May 5th by rotating the crystal 157.5 degrees clockwise (ce — 
157.5 degrees) which corresponds to 9p = 22.5 degrees. The window to detector 
distance (air gap) was 1.4 cm of air. The ORTEC 450 amplifier was used. The PXR 
spectrum is shown in Figure 9. 

Data was obtained on May 8¢ at the same angle for a longer time to improve signal to 
noise. The ORTEC 450 amplifier was replaced by the TENNELEC TC 244 amplifier for 
this and all the subsequent single angle experiments. The amplifier was set for a 12 Us 
peaking time, which gives the best energy resolution. The delay time of the linear gate 
was adjusted based upon observation of the copper foil fluorescence to 30.6 tts. The air 
gap was 1.1 cm. The PXR and calibration spectra are shown in Figures 10 & 11 The 
calibration for May 5‘ is similar to that of May 8" and is not shown. 

Data was also obtained on July 7“ at a Bragg angle of 22.5 degrees. The LINAC was 
tuned to 88.0 MeV and the air gap thickness was 1.0 cm. The amplifier and gating setups 
were the same as those on May 8th. There appeared to be a drift problem, either a drift in 
klystron noise as the machine warmed up, or a detector drift. To reduce the effects of 
this drift the experiment was performed as follows. First a calibration was obtained using 
the sandwich foil starting at 0 hours, then using the carbon crystal, a PXR spectrum was 
obtained for two hours. Due to a setup error the crystal was positioned such that the 
electron beam missed the tin backing, and therefore there wasn't a tin peak for the 
calibration and later for an absolute photons per electron calculation. This spectrum was 
saved and a second calibration was obtained after three hours. This was followed by two 


separate two hour PXR spectra, which were obtained without moving the ladder. This 
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Figure 9. PXR spectrum from the {002} planes of the carbon graphite crystal, obtained 
on May 5 at an angle of 22.5 degrees. This spectrum clearly shows 6 orders. The low 
energy peak is LINAC induced noise. 
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Figure 10. PXR spectrum from the {002} planes of the carbon graphite crystal, obtained 
on May 8 at an angle of 22.5 degrees. This spectrum clearly shows 7 orders. The pulse 
height analyzer discriminator was set at approximately 4 KeV to discriminate out the 
LINAC induced noise. 
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Figure 11. Energy calibration spectrum obtained on May 8. _ This calibration used the 


spectrum of three separate foils which were added together by use of a spreadsheet on 
QPRO. The pulse height analyzer discriminator was set at approximately 4 KeV to 


discriminate out the LINAC induced noise. 
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resulted in three independent PXR spectra at a nominal Bragg angle of 22.5 degrees. The 
crystal was then rotated to a nominal Bragg angle of 23 degrees and a spectrum was 
obtained for two hours, followed by a third calibration after 9 hours from the beginning ot 
this day's data accumulation. 

The calibration information is shown in Table 2. This table shows that the second 
calibration's centroids appear to shift about six channels to the left compared to their 
corresponding peaks in the first calibration, while the third calibration's centroids appear 
within one channel of the corresponding peaks in the first calibration. This information 
shows that it is crucial to have tin or another fluorescent foil in place for calibration while 
taking a PXR spectrum. It also suggests that data be taken in several shorter time periods 
rather then one longer time period. With knowledge of the tin peak's position for each run 
one can then shift the spectra appropriately and then add them together. This procedure 
will markedly improve the spectral resolution. 

The spectra taken at 22.5 degrees were combined in the following manner. The centroid 
of the tin calibration peak in each spectra was determined. Then by shifting the channels of 
the three spectra so that the tin centroids were the same, the three spectra were added 
channel by channel on a QPRO spreadsheet [Ref. 22]. Figures 12-14 show the combined 
PXR spectrum at 22.5 degrees, the spectrum at 23 degrees and the energy calibration 
spectrum respectively. 

On July 8t similar data was obtained, using the new linear gate set at 35 ts wide and 
delayed from the machine start cycle time, To by 7.7 us. The air gap thickness was 1.1 
cm. Additionally the amplifier peaking time was reduced from 12 Us to 8 us. A lower 
peaking time reduces the effect of double counting single events since the effects of double 
counting are proportional to the pulse peaking time, however there is a tradeoff with 
energy resolution. The energy resolution is related to the pulse peaking time. There are 
two contributions. One is from parallel noise (detector leakage current, and all resistors 


which parallel the detector) that is proportion to the square root of the peaking time. The 


a7 


TABLE 2. JULY 7h CALIBRATION DATA, SHOWING PEAK ENERGY, 
ENTROP CHANNEL NUMBER AND PEAK FWHM. 
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Figure 12. PXR spectrum from the {002} planes of carbon graphite, obtained on July 7th 
at an angle of 22.5 degrees. This spectrum was obtained by combining three separate 
Shorter time period PXR spectrum. This spectrum shows six orders with the fifth order 
and the tin calibration peak slightly overlapping, with the tin peak at a higher energy than 
the n=5 peak. 
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Figure 13. PXR spectrum from the {002} planes of carbon graphite, obtained on July 7th 
at an angle of 23. degrees. This spectrum shows six orders with the fifth order and the 
tin calibration peak slightly overlapping, with the tin peak at a higher energy than the n=5 
peak. 
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Figure 14. First of three energy calibrations obtained on July 78. This calibration used 
the sandwich foil containing Titanium, Yttrium and Tin. The Kg and Kp peaks of Y and 


Sn are both visible, while the Ti only shows the Kg peak. 
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second is serial noise (total input capacitance, transconductance (gain of the input device 
etc.) ) that is inversely proportional to the square root of the peaking time. This results in 
an optimum peaking time for SiLi detectors between 8 and 24 us. [Ref. 20] The LINAC 
was tuned to 84.5 MeV. Additionally, the PHA was updated to the PHA-II software 
which has 8192 channels compared to 2048 for PHA-I. The PHA-II digital stabilization 
was utilized to counter the effects of drift. An energy calibration was obtained before and 
after the PXR spectrum. To compare the results obtained with the results from the 2048 
channel analyzer the data was binned using a simple program written in C (see Appendix 
B). Four channels were added together and placed in one bin using the energy value 
obtained from the third of the four bins. Figures 15 & 16 show the binned PXR spectrum 
and the first binned energy calibration spectrum. The resolution is probably improved 
because of collecting data over more channels and then binning, and due to the use of the 


PHA-II digital stabilization option. 


B. SINGLE ANGLE PXR ANALYSIS (FAR REGION) 

In order to compare the data with theory, and from one experiment to the next the 
following analysis was performed. First the net area of each spectral peak was 
determined. This was done two ways. 

The first way used the software supplied with the PCA I TENNELEC PHA. Using the 
PHA computer analysis the peaks were determined by estimating where each peak began 
and ended, and setting these channels as a region of interest. The PHA analyzer 
determines the average background by averaging the start channel with the three previous 
channels and the stop channel with the three following channels and drawing a straight line 
connecting these two average values. All counts below this noise line are treated as noise 
and counts above it as signal. The program determines the gross and net area. Gross area 


is the total number of counts in the region (including noise) and net area 1s only the 
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Figure 15. PXR spectrum from the {002} planes of carbon graphite, obtained on July 8th 
at an angle of 22.5 degrees, using the Canberra (SiLi) detector. This spectrum shows six 
orders with the fifth order and the tin calibration peak slightly overlapping, with the tin 
peak at a higher energy than the n=5 peak. 
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Figure 16. Energy calibration obtained on July 8th. This calibration used the sandwich 
foil and the Canberra (SiLi) detector. This spectrum shows the Kg peak of Ti, Y and Sn, 


and the Kg peaks of Y and Sn. This spectrum was plotted by binning four channels 


together and plotted at the central channel of these four bins. 
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number of counts above the noise line. It also determines the FWHM and the centroid of 


the peak using the following equation: 


>) counts, cnn 


of 
“ ) counts, 
1 


(31) 


where Chn; is the channel number (or equivalent energy) and counts; 1s the corresponding 
number of counts in that channel. 

The second way was to determine the average noise in the regions between the peaks 
using the gross area of 100 channels divided by 100. The background was determined by 
multiplying the number of channels in the peak by the average noise in the region just 
following the peak. The net area is determined by subtracting this background from the 


Gross area. The error in the Net area 1s 


AN o) ATCA ocoss + rater a entre 
— CO T—__ 
N BC Aree CA aun 


gross 


(32) 


since the error in the area is VArea. Tables 3-7 show these values for the various 
experiments. Figure 17 shows a plot of the comparison between the PHA net area and the 
Extended Background net area for the May 8th data. This plot shows that except for the 
first two peaks there 1s very little difference between the net areas as determined by these 
two methods. | 

The raw spectrum area needs to be corrected for attenuation as the photon travels 
between the target and the detector. Also the relative efficiency of the detector must be 
taken into account. The attenuation coefficients were determined using a computer 
program called XCOM [Ref 21] near the centroid energy. Some of the coefficients are 
listed in Table 8, along with the detector efficiency [Ref. 5] and the values of the 


correction factors for Be, Kapton and for air at the two extreme thickness values. The 


final column shows the total correction factor for an air thickness of 1.1 cm. 
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TABLE 3: MOSAIC CARBON CRYSTAL{002} PLANE PEAK DATA TAKEN MAY 
sth FOR 90 MEV ELECTRON BEAM AT BRAGG ANGLE OF Op = 22.5 DEGREES 
FOR PEAKS 1-6. 
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TABLE 4: MOSAIC CARBON CRYSTAL{002} PLANE PEAK DATA TAKEN MAY 
gth FOR 90 MEV ELECTRON BEAM AT BRAGG ANGLE OF 6p = 22.5 DEGREES 


FOR PEAKS 1-4. 
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TABLE 5. SAME AS TABLE 4 FOR PEAKS 5-8. 
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TABLE 6: MOSAIC CARBON CRYSTAL{002} PLANE COMBINED PEAK DATA 
TAKEN JULY 74 FOR 88 MEV ELECTRON BEAM AT BRAGG ANGLE OF 68 = 
22.5 DEGREES FOR PEAKS 1-5. ANALYZED BY COMBINING THE THREE 
INDIVIDUAL 22.5 DEGREE SPECTRUM. 
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TABLE 7. MOSAIC CARBON CRYSTAL{002} PLANE PEAK DATA TAKEN JULY 
7th FOR 88 MEV ELECTRON BEAM AT BRAGG ANGLE OF 6g = 23.0 DEGREES 
FOR PEAKS 1-5. 
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[The net area was determined by use of a peak fitting program with the following parameters: Gaussian 
Peak: 25.197 KeV centroid, 0.548 FWHM with a background of counts = 20.94 - 0.527*Energy, which 1s 
consistent with the balue of 7 used in the Peak channel. 
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Plot of the net area in each peak as a function of the peak number as 


determined by the pulse height analyzer and by using the extended background method for 
4] 


the May 8th data. 


Figure 17. 


TABLE 8 ATTENUATION COEFFICIENTS FOR THE PEAK ENERGIES, 
DENSITY AND THICKNESS FOR BE, DRY AIR AND KAPTON. 
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Using the net area from the PHA analysis and the coefficients listed above the corrected 
counts 1s determined using the following: 

= [ap Tet Hea HOt (33) 
where I is either Net Area or Net Peak counts counted by the (SiL1) detector. 

Two methods of analysis were used to determine the counts at the front face of the 
crystal Ig. The first and simplest used equation 33 applied peak to peak where (UPt)total 
is assumed constant over the entire peak. This is the analysis used in references 4 & 5. 
However, for the low energy peaks p varies over the peak and the values obtained for Io 


vary depending on if you use {1(5 keV) or u for each individual channel. For this reason a 


second method was examined. The second method involved using a least squares fit of 
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the attenuation coefficient. The actual coefficients and the approximations for the three 
materials are shown in Figures 18-20. The formulas give U as a function of energy and 
are applied to the values of energy and the number of counts on a QuattroPro (QPRO) 
spreadsheet. [Ref. 22] The counts are then summed over the peak and a new Io is 
obtained. See Table 9 for the areas and the ratio of the area of the nth order to the area of 
the first order for the raw data as collected by the detector, the corrected data using the 
linear attenuation coefficient at the centroid and the data corrected using the linear 
attenuation coefficient for each channel. Figure 21 plots the raw area, the area corrected 
for attenuation based on the coefficient at the peak centroid and the coefficient for each 
channel using the least square fit values for the May 8th data. 

It is not clear that the use of the channel by channel correction is appropriate since the 
detector response not only broadens the PXR spectrum, there is a folding of the detector 
response into the PXR spectral response. If the detector resolution is small compared to 
the PXR resolution the channel by channel correction method is a better approximation, if 
the PXR resolution is small compared to the detector resolution the correction at the 
centroid is the better approximation, and if the two are comparable which is the case for 
these experiments as seen below the best method would be to deconvolve the PXR 
spectrum from the detector response and then apply a channel by channel correction. 
Since this method 1s extremely complex it is not applied here, and for ease of calculation 
the method of using the attenuation coefficient at the peak centroid will be utilized 
for the remainder of the analysis. In order to determine the best method the PXR 
resolution needs to be calculated. 

1. Detector Resolution at PXR Spectral Energies. 

To determine the detector energy resolution, the calibration peaks were assumed to 


be Gaussian in shape. The natural line widths of the x-ray fluoresence peaks were 
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Figure 18. Plot of the log of the actual linear attenuation coefficient for Be vs. energy in 
KeV, and the formula values using 3 section least square fit. 
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Figure 19. Plot of the log of the actual linear attenuation coefficient for Dry Air vs. energy 
in KeV, and the formula values using 3 section least square fit. 
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Figure 20. Plot of the log of the actual linear attenuation coefficient for Kapton vs. energy 
in KeV, and the formula values using 3 section least square fit. 


obtained from Reference 23, these natural line widths varied from 10 eV to 45 eV. The 
natural line widths were also assumed to be Gaussian distributions. From the measured 


fluorescence peaks the detector resolution or linewidths were determined by adding them 


in quadrature: 


ae ee (34) 


obs det lw? 


where I,,is the measured linewidth, I°,., is the detector resolution and I, is the natural 
fluorescent linewidth. The units:are keV. The uncertainty in I, is assumed to be zero, 
therefore the uncertainty in [,., is the same as the uncertainty in T.,,.. The experimental 


values of IT 


det 


and I,. by less than two tenths of a percent. These results are in 


agreement with those previouly obtained [Ref. 4.] Thus the value of the detector 


resolution to be used will be that of the observed calibration spectrum line widths. 
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TABLE 9. DATA CORRECTED FOR ATTENUATION THROUGH KAPTON 
WINDOW, AIR AND BE DETECTOR WINDOW. 


pe | 1 | 2 | 3 |e 


Area Corrected for attenuation at the peak centroid; ratio Area,/Area 


8 Ma 30124; 1 {31742 ;: 1.05,15575 ; 0.52) 7379 ; 0.24 | 3863 Gai 
rea corrected at each channel using fitted coefficients; ratio Area,/Area 
rsuly22.s° | 248161 [22484 090) 0498 -038| 4025:016| | 
rsuy23’ __|_20843-1 | astosa1fusss7.oss| soeio2a | 
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Figure 21. Plot of the raw area at the detector, the corrected area using the linear 
attenuation coefficient at the centroid of the peak (labeled centroid) and the linear 
attenuation coefficient for each channel based on a least square fit for the May 8th data. 


Table 10 lists the values of the natural x-ray fluorescence linewidths and the observed 


linewidths as determined by the pulse height analyzer for the calibration foils for the May 


Sth May 8" July 7th and July 8' spectra. 


TABLE 10. LIST OF CALIBRATION FOILS AND THE RESPECTIVE X-RAY 
ENERGY PEAK VALUES. INCLUDED IN THE LIST ARE THE NATURE 


LINEWIDTHS OF THE FOILS, AND THE OBSERVED RESOLUTION. 


Po ee 


ws | cu | wos | 003 | oe 

a ee 
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a ee ee 
a ee ee ee 
a ee ee ee 
et vf taf 06 | 0300 
a ee ee ee 















-The linewidth of Y wasn't listed in reference 23. The values was approximated by using the values of Sr 


(56:2 and Zr (0.9 eV). 


3The linewidth of In wasn't listed in reference 23. The value was approximated by using the values of Pd 


(10.2eV) and Ag (10.0 eV) . 
+The observed resolution was taken from the calibration taken after three hours. 


’The linewidth of Sn wasn't listed in reference 23. The value was approximated by using the values of Pd 
& Ag (10 eV). The next element listed was W (43 eV) giving an upper limit which is sull insignificant 


compared the the observed resolution. 
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2. PXR Bandwidth Determination. 


The PXR bandwidth is determined in the same manner, that is 


me. = [,, 


det? 


(35) 
where I’pyp is the bandwith of the PXR spectral line. The units arein keV. Table 11 lists 
the values of the observed resolution for a given energy, the calculated PXR bandwidth 
for the various single angle experiments and the expected value as calculated using 
equation (27). 

Note that the observed linewidths in Table 10. vary with the various experiments. 
These variation are due to the changes made in the experimental procedure. That is the 
May 5th data is much wider due to poorer resolution of the ORTEC 450 amplifier. This 
amplifier was replaced by the TENNELEC TC 244 amplifier following that experiment. 
The July 8th data used an upgraded version of the PHA software which has _ digital 
stabilization features which compensate for drift and therefore reduces the observed 
linewidth. Notice that in spite of these changes in detector resolution, the PXR linewidths 
do not vary significantly, since the changes in the detector resolution are accounted for by 
equation (35). 


TABLE 11. LIST OF OBSERVED RESOLUTION AND CALCULATED PXR 
BANDWIDTH FOR THE ow SINGLE ANGLE EXPERIMENTS. 
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TABLE 11. cont. 


| Date | Energy (keV) | Typ. (keV) | Tpxp (keV) 


degrees 
en nn eas 
| t—“‘zLS 43 | 063 | 0 
TS 

0.87 

=e 

degrees 
———- 
ey os 
CT S052 ee eee 
| rt—“‘zL]SSCUCL 2G | OO ee 


Once the resolution and the net area is determined at the front face of the crystal, 












these values can be compared to the theoretical values. Since, for most of these 


experiments, there was no reference tin peak from which to compute the LINAC current, 


the analysis is accomplished by comparing the experimental ratios of the intensity (area) of 
a given order to the intensity of the first order with that calculated from theory. The single 
angle experiments were all obtained in the far region, therefore the theoretical ratios are 
determined by integrating equation (26) over the field of view, and the bandwidth 1s based 
on A@, and not the FWHM of Jo. The experimental ratios are listed in Table 9 above. As 
described in Chapter II , section C, the calculations must be performed over the limited 
field of view. 

Figure 2 in Chapter II section C illustrates how the aperture in the far case 1s less 
than the FWHM of J>_ The theoretical ratios are calculated by integrating equation (26) 
over the field of view for each order then taking the ratio of the number for a given order 


to the first order. Since the effects of mosaic spread and multiple scattering aren't well 


6The observed resolution was taken from the calibration taken after three hours. 
Estimated due to overlap of fifth order peak with the tin peak in the PXR spectrum. 
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understood, the theoretical values with and without these effects are included for 
comparison. These ratios are shown in Table 12 along with the experimental ratios. The 
experimental ratios are plotted as a function of order in Figure 22. The experimental 
values are compared with the theory by dividing the experimental value by the theoretical 


value and is also shown in Table 12. 


TABLE 12. COMPARISION OF EXPERIMENTAL AREA RATIOS, I[(n)/I(1), WITH 
THEORY WITH AND WITHOUT THE EFFECTS OF MULTIPLE SCATTERING 
AND MOSAIC SPREAD. 


| Experimental Data 
















Theory with 
scattering and 
mosaic spread 


| spread | 
Date | nj] Net Area Ratio he (| Ratio | Exp 
! (counts) theory theory } 
| hs a 
| aes [Poo ia 
| ees 
aes 
| 5) 
ws 02 00mm 1 
ee Lit ies 
ee | 
| aE 
| ll 
| 
[009 | 9. 
aE [00s | 68 
aE 
an a 
a. | 3 oO? es 


Theory without 
scattering or mosaic 
























KES 
ee 4 
| a 
8 July | 1 
me | 2 
| ED Jato oa AO) 
| Es 
| ES 0.005 | _80_ 









2) 


SS: 


PLP LPL LMM a: 
www 


ssssssssssscsssssss258% 
ae, 


a 


eee eee e ee = 
Bee seu BS 
Sao S22 SES SSE S2SS8 SSS 8S S288288 828282 8828288282822 8:' 


Ratio Ieny/IO) 
SASS 


a 


SSNS: 


= 


“lita 


eH 7 
My ; oo 
a Ve 


Hm ONE 
| Nin Gs 


? \' 
y | \' 
aN 
iN 
WN 
NN 
AN 
a \ 
LAN 
TAY 
Ww 
: N N“ 








May 5 May 8 July 7 22.5 degree 


Eee July? 23 degrees PXd July 8 






Figure 22. Plot of Area ratio of nth order peak to first order peak vs order for the four 
experiments conducted near the Bragg angle of 22.5 degreees. 


Notice that the experimental ratios are consistantly larger than theory, and the inclusion of 


mosaic spread and multiple scattering makes the discrepancy even larger, thus showing 


a2 


that equation (21) may be preferred over equation (20). This indicates that the theory is 
inadequate to explain the observed results. It also indicates that it is important to look at 
the convolution of the electron beam divergence with the multiple scattering effect for an 
even better estimation of the effect of multiple scattering. 

3. Calculations of Absolute PXR Yield (Photons/Electron) 

In addition to comparing the intensity ratios, for the 7 and 8 July experiments the 
number of photons per electron can be calculated. The Linac current can be determined 


by using the yield of the Tin fluorescent peak using the following equation 


Afedas 
[ lime | AW) 


Luinac = —— ey (36) 
Oia 2% 4 JOON, Butea) 

where Ij -NAC 1S the calculated beam current, Areay;, 1s the Area of the tin peak at the 
detector (from the pulse height analyzer), time is the collection time of the PXR spectrum, 
AW is the atomic weight, o is the electron interaction cross section (1.0 e-26 m”)[Ref. 
24], Qget is the solid angle subtended by the detector (2x10-* sr.), p is the density of the 
target (2.260 g/cm?), Ng is Avogadro's number, fgey is the de-excitation transition 
probability (0.712)[Ref. 25], t is the thickness of the target (2.98x107>m for 22.5 degrees), 
€ 1s the relative detector efficiency at the tin peak energy (0.9 for these experiments) and a 
is the total attenuation factor of photons (through the effective thickness of tin, the 
effective thickness of the mosaic crystal, the kapton window, the air gap and the beryllium 
window of the detector (0.763 for 22.5 degrees. )[Ref. 21]). 

The electron interaction cross sections [Ref. 24] were obtained for 85 MeV 
electrons. The transition probabilities [Ref. 25] include both the K-L11 and the K-L111 
de-excitation probabilities. The attenuation coefficients were obtained using XCOM at the 
known de-excitation energy of Tin. The details of these calculations may be found in 


reference 17. 
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Using the LINAC current in electrons per second, the intensity of the PXR in counts 
as listed in Table 9 and the data collection time, the number of photons can be calculated 
as; 

(PXR Area, 


hoe 


8 (37) 
LINAC Current 


These values are shown tn Table 13 along with the calculated currents. 


TABLE 13. COMPARISON OF THE EXPERIMENTAL NUMBER OF PHOTONS 
PER ELECTRON TO THEORY WITH AND WITHOUT THE EFFECTS OF MOSAIC 
SPREAD AND MULTIPLE SCATTERING. 
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The number of photons per electron are plotted vs order for these three dates in 
Figure 23. Note that the data 1s consistently larger than the theory, as seen with the ratios. 
Figure 24 is a plot of the number of photons per electron for experimental data, and the 
theory with and without the effects of multiple scattering and mosaic spread. With the 
actual number of photons per electron there is more information than with the ratios alone. 
First it is clear that the data is much larger than expected. Second, the difference between 
the data and the theory gets larger with order, which explains why the error in the ratios 
increased with order. Because of overlap, a peak fitting routine was used to obtain the 


area under the fifth order PXR and tin peaks. [Ref. 26] 
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Figure 23. Plot of the absolute number of photons per electron vs order for the three 
experiments near the Bragg angle of 22.5 degrees with a tin peak used to determine the 
LINAC current. 
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Figure 24. Plot of number of photons per electron vs angle for the 7 July 22.5 degree data 
along with the theory both with and without the effects of multiple scattering and mosaic 
spread. 
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C. PXR DATA COLLECTION (DISPLACEMENT FROM THE BRAGG ANGLE) 
Two sets of "rocking data" were obtained. For these experiments rocking data refers to 
a fixed detector angle of 45 degrees with respect to the beam line, and the crystal angle is 


rocked about the Bragg angle of 22.5 degrees. The first was obtained on June 5, 1992. 





Data was obtained using the ORTEC SiL1 detector and the TENNELEC amplifier. This is 





the far region. The linear gate was set at S Us wide with a delay time of 30.6 Us from To 
(the LINAC start sequence time). The air gap thickness was 1.1 cm. Brief spectra were 
obtained starting at 22.5 degrees for 200 seconds. Then the edges were determined by 
obtaining spectra at 26.5 and 18.5 degrees, both of which showed no PXR peaks. To 
obtain an estimate of the angle giving the largest count rate three more angles were set for 
200 seconds each. The results are shown in Table 14, where "counts" is the net number of 
counts in the peak channel of the n=1 peak. The count rate appeared to increase upon 
rotation from 22.5 degrees to lower angles with the maximum near 22.1 degrees. Angles 
were selected assuming 22.1 degrees to be the angle which produced the maximum count 
rate. (This doesn't agree with theory where 22.5 degrees should be the maximum, but due 


to uncertainties in the crystal onentation this was assumed within error). 


TABLE 14. NUMBER OF COUNTS VS. ANGLE FOR 200 SECOND SPECTRUM 
OBTAINED ON JUNE 5, 1992 (FAR REGION). 


ass | 
pws | 
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The crystal angle was varied starting at 19.1 degrees, and increased by 1.0 degree. 
When the angle was within | degree of the assumed maximum the angle was increased by 
0.5 degrees. The spectra were obtained for 30 minutes each. As seen in Table 15 the 
number of counts in the peak channel in the n=1 peak continued to increase beyond 22.1 
degrees. The angle was increased in 0.5 degree steps until the number of counts in the 
peak channel were comparable to the number at 21.1 degrees, after which 1.0 degree steps 
were taken until very little n=l peak showed. These uncorrected PXR spectra are shown 
in Figures 25 & 26. Full scale spectra can be found in Appendix A. This data shows that 
the maximum count rate 1s obtained near 22.5 degrees as expected. 

On July 22, 1992, rocking data was obtained using the Canberra (SiL1) detector and the 
TENNELEC amplifier. This is the near region. The LINAC was tuned to 95.9 MeV. 
The air gap was 1.1 cm. Due to a set up error the data collection time was not fixed and 
therefore varied at different angles but was approximately 45 minutes. A longer time was 
used than on June 5" since the Tin peak wasn't as well defined. The angles where taken 
to be the same as those used on June 5th for comparison purposes with the exception that 
no spectrum was collected at 19.1 degrees. These uncorrected PXR spectra are shown in 


Figures 27 & 28. Full scale spectra can be found in Appendix A. 


D. PXR ANALYSIS (DISPLACEMENT FROM THE BRAGG ANGLE) 

The analysis of the data involved three issues. The first is the value of the PXR energy 
versus angle, the second 1s the intensities versus angle and the last is the number of 
photons per electron versus angle. 

To determine the PXR energy versus angle the original spectra were analyzed using the 
pulse height analyzer to determine the centroid of each peak, including the centroid of the 
tin peak. The value of the first order centroid was plotted vs. angle, along with the 
theoretical value. There appeared to be a discrepancy between the theory and the data. 


Next the tin peak centroid was plotted vs. angle. Note that theoretically the tin peak 
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TABLE 15. NUMBER OF COUNTS VS. ANGLE FOR 200 SECOND SPECTR 
OBTAINED ON JURY 2251732 AR REGION), 
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Figure 25. PXR spectra from the {002} planes of carbon graphite for the "far region", 
obtained on June Sth at angles 19.1-22.6 degrees. 
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Figure 26. PXR spectra from the {002} planes of carbon graphite for the "far region", 
obtained on June Sth at angles 22.6-25.6. 
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Figure 27. PXR spectra from the {002} planes of carbon graphite for the "near region", 
obtained July 22nd at angles 20.1-22.6. 
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Figure 28. PXR spectra from the {002} planes of carbon graphite for the "near region", 
obtaines July 22nd at angles 22 6-25.6. 
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position in the spectrum should not depend on angle. The fact that the tin peak's position 
was not constant with angle showed that there was some shifting in energy with time. 
Taking the difference between the experimental tin peak centroid and the known tin de- 
excitation energy (25.1955KeV) and subtracting this difference from the PXR centroid, 
this corrected PXR centroid was compared with theory and showed much better 
agreement. These plots are shown in Figures 29 & 30 for the June sth and July 22nd 
experiments respectively. The scale on the left refers to the PXR energies. The original 
PXR centroids are plotted as "pxr centroid", the PXR centroid minus the difference in the 
true tin and known tin is plotted as "corrected cent”. The experimental tin peak is shown 
using the scale on the right, with the known tin energy plotted as a dotted line. Note that 
the same increments were used on the two scales for comparison purposes. 

This also explains the difference between the expected values of E(n)/E(1) equal to n 
and the actual values as listed in Tables 2-6. This is further illustrated in Table 16. The 
(*) indicates corrected values based on the shift in the Tin peak centroid energy from the 


known de-excitation energy. This shift 1s listed as "Tin Shift”. 


TABLE 16. COMPARISON OF RAW ENERGY RATIOS TO THE RATIOS OF THE 
SHIFTED ENERGIES FOR THE 5 JUNE ROCKING DATA. 
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Figure 29. Plot of PXR first order energy orginal data, corrected data and theory vs. 
angle, tin centroid and theory for the June 5t4 experiment (far region). 
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Figure 30. Plot of PXR first order energy onginal data, corrected data and theory vs. 
angle, tin centroid and theory for the July 22"4 experiment (near region). 
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These values were computed by finding the energy centroid based on equations (23) and 


(24) within the limits of the field of view. That is: 


2 SES (37) 


where the integrals are computed over the field of view A@.. Figures 31-33 illustrate the 
effect of changing angles on the field of view. Figure 31 illustrates A@, for the far case 
near the Bragg angle, while Figure 32 is approximately one degree displaced from the 
Bragg angle for the far case. Notice that the effect of changing the crystal angle is to shift 
the field of view across the function Jo. The dector sees a slice of J> and that slice shifts 
across J> as the crystal angle changes from the Bragg angle. Figure 33 illustrates the field 
of view for the near case. Notice that the field of view is larger than the FWHM of J> for 
this case. 

Comparison of Figures 29 & 30 show that for the far region the detail of the function Jo 
is not important in the energy centroid calculations, however for the near region the detail 
of Jo has a definite effect on the calculations. The dip in Figure 30 near 23 degrees shows 
the effect of the dip of Jo. There is a larger discrepancy between the theory and the 
corrected data for the July 22nd experiment. This may be due to the small tin peak area 
and therefore a larger error in determining the tin peak centroid. Also, since the details of 
J> have a larger effect, an error in the angle determination has a larger effect, and shifting 
the data by 1/4 degree would show better agreement with the corrected data and theory. 
Figures 34 and 35 shows a plot of the corrected values of the energies for the first three 
orders versus angle, along with the theoretical values for the June sth and July 22nd 
experiments respectively. Note that on this scale the difference between theory and 
experiment is very small, although there is a larger difference at the smaller angles for the 


higher orders. 
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Figure 31. Illustration J> and the field of view for the far region near the Bragg angle of 
22.5 degrees. 
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Figure 32. Illustration of J> and the field of view for the far region approximately one 
degree displaced from the Bragg angle i.e. 23.5 degrees. 
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Figure 33. Illustration of J> and the field of view for the near region near the Bragg angle 
of 22.5 degrees. 
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Figure 34. Plot of PXR corrected data and theory vs. angle for the first three orders for 
the June 5th experiment (far region). 
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Figure 35. Plot of PXR corrected data and theory vs. angle for the first three orders for 
the July 2294 experiment (near region). 
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To compare the intensities at each angle the spectra obtained on June sth and July 22nd 
at various angles needed to be normalized . Since the LINAC current was adjusted to 
maintain a machine pulse to PHA count rate ratio of about 6:1, and since angles closer to 
the Bragg angle of 22.5 degrees produce a larger number of x-rays, the current at these 
angles was lower. 

For both dates, one spectrum was used to "strip" all the other spectra by shifting and 
scaling the spectrum based on the tin peaks. For the June Sth data the 19.1 degree 
spectrum was used and for the July 22nd data the 20.1 spectrum was used. The scaling 
factor was determined by equation 33. 


Tin - Area Ref Ti 
Scale Factor = paca 1 ime Caer 110. (37) 
Area Tin 


where Area Tin is the actual tin area of the spectrum being scaled and Area Ref Tin is the 
corresponding area for the equivalent channels of the 19.1 or 20.1 degree spectrum tin 
peak. To account for the shift in tin energy, the spectrum was shifted and the reference 
was subtracted until the final spectrum showed an approximately zero net area. 

After the spectrum were scaled and shifted, the June 5th data showed an average 
background of zero. For the July 22"¢ data there remained small amounts of LINAC 
induced noise, especially at the lower energies. Using QPRO spreadsheets and the XCOM 
program the counts were corrected, channel by channel for attenuation through the Kapton 
window, the air gap and the Be window of the detector. The intensity (area) at the front 
face of the crystal was determined by summing over the channels of each peak. The areas 
are listed in Tables 16 and 17 for the far region (June Sth experiment) and the near region 
(July 22nd experiment) respectively. The scaled, shifted and corrected spectra for both 
the June 5th and July 224 experiments can be found in Appendix A. 

Using the scaled spectra, the intensity can be plotted as a function of angle for the 
various orders. Figures 36 and 37 show the scaled area intensities vs. angle for the June 


5th and July 2204 experiments respectively. Note that the second order peak is larger than 
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TABLE 16. SCALED AREA INSTENSITIES FOR THE JUNE 57TH EXPERIMERS 
(FAR REGION). 





TABLE 17. SCALED AREA INTENSITIES FOR THE JULY 22ND EXPERIMENT 
(NEAR REGION). 


panos | on [| 2 | 3 | «4 | os 
pms | ses | soo | so fe 
pas | ss | sm} ass} 
-so1 | i950 | aaas | sa90 | 2550 | 900 
ere re er er 
pss | izes | oso | asso | sso | 
2s | so | 240 >| sel 
| on | 26s | is | |e 
| rs | noo | mo | - | .~ | 
| ass | em | = | es ee] 









74 


Intesity (counts) 





‘ See 
-% " ae . = 


occa 
21 215 22 25 23 235 24 245 25 255 
Angle in degrees 





=m Peak 1 —* Peak 2 -*- Peak 3 


"= Peak 4 == Peak 5 -ake- Peak 6 





Figure 36. Plot of scaled area intensity for all orders versus angle for the June Sth 
experiment. 
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Figure 37. Plot of scaled area intensities for all orders versus angle for the July 22nd 
experiments. 
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the first order peak of the June 5‘ data at both 22.5 and 23 degrees, but it is smaller than 
the first order peak for all the angles of the July 22nd experiments . Also note that the plot 
as a function of angle for a given order is broader for the July 2294 experiment than for the 
June 5th experiment i.e., the PXR intensity falls off faster as a function of angle for the far 
region than for the near region. 

The next comparison of the data with theory involved analysis using ratios. For the first 
order only, the intensities were divided by the first order intensity at 22.6 degrees (the 
angle with maximum intensity). Although the theoretical maximum is at 22.5 degrees, to 
compare with the experimental data, the number of photons per electron for each angle 
was divided by the theoretical value at 23 degrees. Figures 38 & 39 show these 
comparisons for the June 54 and July 224 experiments respectively. Note that the June 
5th data appears to be shifted by approximately 0.3 degrees, as shown by the dotted line. 
Also Note that the July 22nd data appears to be shifted in the other direction (compared to 
the June Sth data) by about 0.1 degrees. 

These ratios compared one angle's first order intensity to that of the first order intensity 
at 22.6 degrees. A second comparison involves the ratio of any given order intensity to 
the first order intensity for a single angle. Tables 18 and 19 show these theoretical ratios 
with the experimental data in brackets for the far and near case respectively along with the 
experimental ratios in brackets. Figures 40 and 41 plot the experimental and theoretical 
values for I(2)/I(1) and for I(3)/I(1) for both the far and near region respectively. Note 
that there 1s little vanation of the theoretical ratios with angle and that there is little 
difference between the theoretical values for the far and near case. The fact that the 
experimental values vary with angle with a maximum at Op shows that searching for the 
maximum in this ratio can be used as a means to position the crystal at Op if the LINAC 
current isn't monitored. 

Lastly, the number of photons per electron for the data at each angle can be compared 


to the theory. The method of analysis is the same as that used in the single angle 
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Figure 38. Plot of scaled area intensity divided by the area of 22.6 degrees for the first 
order versus angle, the theoretical ratio versus angle and the shifted theoretical ratio for 
the June Sth experiment (far region). 


78 


0,8 


0,6 


Ratio 


0.4 


= 





0 : 
eeeecWm@rel.| lio eel ee 0 ed edo ¢4.1 246 25.1 256 


Angle 





Figure 39. Plot of scaled area intensity divided by the area of 22.6 degrees for the first 
order versus angle and the theoretical ratio versus angle for the July 224 experiments 
(near region). 
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TABLE 18. RATIO OF INTENSITY TO THE INTENSITY OF THE FIRST ORDER 
FOR VARIOUS ANGLES FOR THE JUNE 5th EXPERIMENT. 


ANGLE 1(2)1(1); [DATA] | 1(3)/I(1); [DATA] | 1(4)10.); [DATA 






20. 0.184 0.034 0.009 
211 0.185: [0.51 0.034 0.009 
216 0.188: [0.44 O54 (a2 0.009 
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TABLE 19. RATIO OF INTENSITY TO THE INTENSITY OF THE FIRST ORDER 
FOR VARIOUS ANGLES FOR THE JULY 2294 EXPERIMENT. 


ANGLE 1(2)/1(1); [DATA] | 1(3)1(1); [DATA] | 1(4)/10.); [DATA 







nor | ts 0.034 0.009 
Snr | ose joasy | oosagory | 0.009 
ae 0.198: [0.69 0.036; [0.24 0.010; [0.09 
er | oaraposoy | o.osojossy | oor porm 
ain 0.212; [0.82 0.039: [0.36 0.011; [0.16 








Zou OZ Ons 0.040; [0.37 CLO One 
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Figure 40.Plot of the experimental and theoretical ratio of the nth order intensities to the 
first order (I(n)/I(1)) for the n=2 and n=3 peaks for the June Sth experiment (far region). 
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Figure 41. Plot of the experimental and theoretical ratio of the nth order intensities to the 
first order (I(n)/I(1)) for the n=2 and n=3 peaks for the July 22nd experiment (near 
region). 
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experiments. The current was determined by using the area of the tin peak in the 19.1 and 
21.1 degrees spectra for the June 58 and July 22" experiments respectively. These angle 
were used since the scaling is based on them, and therefore the current will be normalized 
to this current. Table 20 and 21] show the experimental values, the theoretical values and 
the ratio of experiment to theory for the first order peaks at the various angles for the June 
5th and July 22"4 experiments respectively. The theoretical values are calculated using the 
theory presented in Chapter I. The large variation between the far region and the near 
region are due to the larger solid angle subtended by the detector in the near region. A 
larger solid angle corresponds to a larger field of view. These values were computed 
using computer integration of equation (23). An example of the computer variables is 
presented in Appendix B. Again, the theoretical values are much lower than the data in 
both the far and the near case. All of the analysis consistently shows that the data is 


clearly much larger than expected. 


TABLE 20. NUMBER OF PHOTONS PER ELECTRON FOR THE DATA, THEORY 
WITHOUT MOSAIC SPREAD OR MULTIPLE SCATTERING AND THEORY WITH 
MOSAIC SPREAD AND MULTIPLE SCATTERING, FOR THE FIRST ORDER 
PEAK VERSUS ANGLE FOR THE JUNE 5th EXPERIMENT (I = 3.39E-19 AMPS). 


ANGLE | DATA | THEORY DATA THEORY DATA 
WITHOUT | THEORY WITH THEORY 
SCATTERING SCATTERING 
pa | ever | saseo | 10 | sro | is 


ZN 2.16e-7 2.13e-8 10 1.30e-8 l 











22.1 4.64e-7 7.99e-8 3.43e-8 
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TABLE 21. NUMBER OF PHOTONS PER ELECTRON FOR THE DATA, THEORY 
WITHOUT MOSAIC SPREAD OR MULTIPLE SCATTERING AND THEORY WITH 
MOSAIC SPREAD AND MULTIPLE SCATTERING, FOR THE FIRST ORDER 
PEAK VERSUS os FOR THE JULY 2294 EXPERIMENT (I = 4.17X10-!5 
AMPS). Noor \ 


ANGLE DATA THEORY DATA THEORY DATA 
WITHOUT THEORY WITH THEORY 
a SCATTERING 


le ee 6 3.67e-8 













2.44e-8 









].32e-7 4.17e-8 
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Lastly, since the experiment in the near region was done as "rocking data" there wasn't 






a bandwidth resolution calculation computed for the near region in section IV. B.1. Table 
22. lists the values of the natural linewidths of the calibration foils and the observed 
resolution for the July 22nd experiment. Table 23. lists the values of the observed 
resolution for a given energy, the calculated PXR bandwidth and the expected value as 


calculated for the near region. 
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TABLE 22. LIST OF CALIBRATION FOILS AND THE RESPECTIVE X-RAY 
ENERGY PEAK VALUES. INCLUDED IN THE LIST ARE THE NATURAL 
LINEWIDTHS OF THE FOILS, AND THE OBSERVED RESOLUTION. 


ee ee ee 


TABLE 23. LIST OF OBSERVED RESOLUTION AND CALCULATED PXR 
BANDWIDTH FOR THE NEAR REGION FROM THE 22.1 DEGREE PXR 
SPECTRUM OBTAINED JULY 22ND (NEAR REGION). 


Pa (KeV) | Tag (KEV) 


S02 0.450 0.420 


| 8 | 0500 0384 0387 
wen | ozs | ts 















8Unable to determine due to poor statistics in peak. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

Enhanced higher order PXR spectra from the {002} planes of carbon have been 
observed. The ratio of the higher order peaks to the first order peaks (I(n)/I(1)) were 
examined. These ratios were found to be much larger than expected, and the difference 
increases with increasing order. The absolute intensity was determined based on the 
LINAC current calculations from the fluorescent x-ray yield of a 0.0275 mm thick tin foil 
placed behind the graphite crystal. The absolute intensity was found to exceed the 
theoretical values by a minimum factor of approximately 10 when the effect of 
multiple scattering of the electron beam and the effect of the mosaic spread of the 
crystal are not considered. This factor increased to more than 3000 when these 
effects are included. 

This is the first time that the theoretical effects of the limited aperture were examined. 
The effect of varying the crystal angle with a fixed detector angle (with respect to the 
electron beam) was examined. The variations with angle of I(n,8)/I(n,22.6°) are consistent 
with theory when the data and theory are normalized to a given angle (in this case 22.6 
degrees.) However, the absolute yield still greatly exceeds the theoretical values. 
That is, for both the far and the near regions 1(1,6)44ta/I(1,@)theory>> 1 1(n,6)/1(1,8) is 
not constant with respect to angle as is the theory. The maximum in this ratio occurs at 9 
g. If current is not monitored, look for the maximum in this ratio to insure the crystal 
angle is 8p. If current is monitored look for the maximum intensity in the first order peak, 
I(1), to insure the crystal angle is Op. 

Additionally, the PXR spectral resolution was calculated for the far and near region. 
The detector subtended a solid angle of 2x10-4 sr. for the far region and 2.37x10-3 sr. for 
the near region. For the far region the PXR spectral resolution was in agreement with 
theory for the lower orders, but was narrower than expected for the higher orders. For 
the near region the PXR spectral resolution was in agreement with theory, such that the 


resolution was wider for a given energy in the near region than in the far region. 
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Further studies need to be done to further examine the effects of multiple scattering and 
mosaic spread. By comparing a thick pure crystal such as Si with other mosaic crystals 
such as carbon, the effect of mosaic spread can be clarified further. 

This study has also shown that although alignment is not extremely critical it is 
important to take great care in focus and steering of the electron beam. Additionally, the 
effects of detector drift, or shifts due to other causes should be carefully examined. Many 
short data runs should be collected. These can then be added on a spreadsheet after 
shifting the known calibration peaks centroids so that they overlap, thus minimizing these 
effects. 

It 1s extremely important for calibration purposes to have an x-ray fluorescent foil 
behind the crystal. This not only greatly improved the analysis, it allowed for the 
calculation of the LINAC current and the subsequent calculation of the absolute number 
of photons/electron. Care should be taken in selecting a fluorescent foil whose x-ray de- 
excitation energy falls between two PXR spectral peaks. Tin worked fairly well for 
carbon at shallower angles, but overlapped the 5th order peak near the Bragg angle, 
complicating the analysis. 

PXR from a thick mosaic carbon graphite crystal has been shown to be a spectrally 
bright source of narrow bandwidth at energies equal to approximately 5 KeV and integer 
multiples thereof. Further study with other crystals could provide additional sources of 


PXR at other energies as well. 
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APPENDIX A. FIGURES OF ROCKING DATA 
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Figure 42. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 19.1 degrees. This spectrum shows negligible PXR peaks and a tin peak at 
about 25 keV. 
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Figure 43. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 20.1 degrees. This spectrum shows a small first order PXR peak and a tin 
peak at about 25 keV. 
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Figure 44. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 21.1 degrees. This spectrum shows the first 2 orders of PXR peaks and a 
tin peak at about 25 keV. 
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Figure 45. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 21.6 degrees. This spectrum shows three orders of PXR peaks and a tin 
peak at about 25 keV. 
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Figure 46. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 22.1 degrees. This spectrum shows possibly four orders of PXR peaks and 
a tin peak at about 25 keV. 
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Figure 47. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 22.6 degrees. This spectrum shows six orders of PXR peaks. The fifth 
order peak and the Tin calibration peak are overlapping, with the Tin peak at a higher 
energy than the n=5 peak. 
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Figure 48. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 23.1 degrees. This spectrum shows five orders of PXR peaks. The fifth 
order peak and the Tin calibration peak are overlapping even more than at 22.6 degrees, 
with the Tin peak at a higher energy than the n=5 peak. 
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Figure 49. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 23.6 degrees. This spectrum shows possibly four orders of PXR peaks and 


a tin peak at about 25 keV. 
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Figure 50. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 24.1 degrees. This spectrum shows three orders of PXR peaks and a tin 
peak at about 25 keV. 
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Figure 51. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 24.6 degrees. This spectrum shows two orders of PXR peaks and a tin 
peak at about 25 keV. 
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Figure 52. PXR spectrum from the {002} planes of carbon graphite, obtained on 5 June 
at an angle of 25.6 degrees. This spectrum shows possibly two orders of PXR peaks and 
a tin peak at about 25 keV. 
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Figure 53. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 20.1 degrees. This spectrum shows two orders of PXR peaks and the Tin 
Ko and Kg peaks 
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Figure 54. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 21.1 degrees. This spectrum shows the first and second order PXR peak 
and the Tin Kg and Kg peaks. 
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Figure 55. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 21.6 degrees. This spectrum shows the first 3 orders of PXR peaks and the 
Tin Kg and Kg peaks. 
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Figure 56. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 22.1 degrees. This spectrum shows four orders of PXR peaks and the Tin 


Kg and Kg peaks. 
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Figure 57. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 22.6 degrees. This spectrum shows five orders of PXR peaks with the Tin 
peak and n=5 peaks overlapping. 
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Figure 58. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 23.1 degrees. This spectrum shows five orders of PXR peaks. The fifth 
order peak and the Tin calibration peak are overlapping. 
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Figure 59. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 23.6 degrees. This spectrum shows possible five orders of PXR peaks. The 
fifth order peak and the Tin calibration peak are overlapping . 
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Figure 60. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 24.1 degrees. This spectrum shows possibly four orders of PXR peaks and 
a tin peak at about 25 keV. 
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Figure 61. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 24.6 degrees. This spectrum shows three orders of PXR peaks and a tin 
peak at about 25 keV. 
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Figure 62. PXR spectrum from the {002} planes of carbon graphite, obtained on 22 July 
at an angle of 25.6 degrees. This spectrum shows two orders of PXR peaks and the Tin 


Ko and Kg peaks. 
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APPENDIX B 


A. PROGRAM IN C. 

This program was written in C to convert the ASCII output from the PHA-II software 
to a file readable by QPRO. It binned every four channels into one bin and assigned the 
value of the third bin's energy to that bin. This program allowed the 8192 channels to be 
displayed in 2048 channels for comparison with the spectra obtained using the 2048 
channel PHA-I software. To prepare the file, an editor was used to eliminate the commas, 
remove the heading comments and add "2045" to the top of the file. The output was in 
the form of "New channel number, Energy, counts” (comma's included). The parameter 
for the conversion from channel to energy must be changed from experiment to 
experiment depending upon the calibration for that data. See Chapter III section C 


equations (28)-(30). 


#include <stdio.h> 

#include <math.h> 

main() {int ea, eb, ec, ed, 1, N, counts: 
float ca, cb, cc, cd, E; 

int rola, roib, roic, roid; 

scanf("%d" ,&n); 

E=0.0; 

counts=0.0; ca=0; cb=0:; cc=0; cd=0; 
for (i=O0:i<N:;++1) { scanf("Yod%f%d" &ea,&ca,&roia); 
scanf("%od%f%d" &eb,&cb,&roib): 
scanf("Yod%f%od" &ec, &cc,&roic): 
scanf("Yod%f%od" &ed &cd,&roid); 


counts=(int)(cat+cb+cc+cd); 
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E=0.00472*(1*4+2)-1.04; 
printf("%d",1); printf(","); 
printf("%g",E); printf(","); 
printf("%od" counts); printf("\n"); } } 


B. COMPUTER INTEGRATION PROGRAM EXAMPLE 

This program was used to determine the absolute intensity and the energy bandwidth 
for the PXR spectra. The various parameters changed based upon the experimental 
conditions. This example uses the variable for the far region at a crystal angle of 22.5 
degrees. It includes the effects of multiple scattering and mosaic spread. 
ne | 
enbeam=90.0 
flvec=(2.01,1.34,0.0946,0.624,0.402,0.268) 
fl=flvec(n] 
z=6 
an=12 
thetamm=0.4 
thetabb=22.5 
dang=6.708 
Lavec= (0.0213,0.181,0.573,1.104,1.586,1.939) 
La=Lavec([n] 
dthetax=15/1002.19 
dthetay=15/1002.19 
dnm=dang/10 
gam=(.511+enbeam)/.511 
thetam=(pi/180)*thetamm 


thetab=(pi/180)*thetabb 
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L=.139/sin(thetab) 

num=La“%2-(I*la+La‘%3)*exp(-L/La) 

den=La*(1-exp(-L/La) 

Tave=num/den 

k=0.157*z*(z+1)*Tave/An*enbeam”2) 

Thetassqr=k*Ln(1.13 *z*(4/3 )an*(-1)* Tave) 
energy=n*2*pi*(197.32)/(dnm*sin(thetab)) 

lambcm=0.000124/energy 

lambang=lambcm* 100000000 

x0=78 192*6/(energy)*2 

x 10sqr=(1.039e-12)* (lambang*4) *f1%2 
thetap=sqrt((1/gam“2)+xo+thetassqr+thetam’2) 

alphay=dthetay/(2 *thetap) 

alphax=dthetax/(2*thetap) 

cs2=(cos(2*thetab))*2 

IN@)H=CEEX 2 “(1 -€s2))/ Cx 2k 

Q(x)=atan(alphay/(1+x%2)*.5) 

j2a(x)=x%2 *cs2*((alphay/(1+x‘%2))* 1/(alphay*2+ 1+x%2)+Q(x)/(1+x%2)*1.5) 
j2(x)=G2a(x)+ Q(x )/(1+x2)*.5-alphay/(alphay*2+ 1+x%2))) 
factor=(1-exp(-L/La)) 

coeff=x 1 Osqr*La/(4*sin(thetab)*cos(thetab)lambcm*thetap *energy* 137)*factor 
dNdE 1(x)=pi* coeff*j1(x) 

dNdE2(x)=2*coeff*j2(x) 

xd NdE2(x)=x*dNdE2(x) 

x 1=((pt/180)*(45-2*thetabb)+. Sdthetax)/thetap 
x2=((p/180)*(45-2*thetabb)-. Sdthetax)/thetap 

xcent=(in((1e-5 *xdNdE2(x)),x=x1 to x2))/(in((1le-5 *dNdE2(x)),x=x1 to x2)) 
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ecent=xcent *energy*thetap/tan(thetab)+energy 
area=in(dNdE2(x), x = xl to x2) 


resolution =(x1-x2)*energy*thetap/tan(thetab). 
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